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^jbstract 

An  experimental  determination  of  the  interframe  buckling  collapse  of  a  ring  stiffened  cylin¬ 
der  subjected  to  hydrostatic  pressure  is  presented.  This  test  served  two  purposes,  the  first  to 
experimentally  verify  finite  element  predictions  of  the  buckling  collapse  behaviour,  the  second 
to  evaluate  the  possibilities  of  using  an  existing  high  pressure  facility  for  conducting  model  tests 
of  more  complex  pressure  vessel  models.  The  finite  element  and  experimental  results  for  the 
collapse  load  and  mode  were  in  good  agreement.  A  method  of  determining  the  elastic  buck¬ 
ling  collapse  load  by  measuring  the  decreasing  frequencies  of  vibration  modes  corresponding  to 
buckling  collapse  modes  has  also  been  investigated  with  encouraging  results.  The  high  pres¬ 
sure  tank  served  the  purpose  of  conducting  the  experiment,  although  the  inability  to  physically 
observe  the  test  did  cause  some  difficulties  in  determining  the  actual  collapse  load  , 

^  0 


Resume 

L’dtude  cxpdrimentale  porte  sur  le  flambage  entre  bads  d'un  cylindre  renforc6  de 
raidisseurs  circulates  et  soumis  k  une  pression  hydrostatique.  L'essai  avait  pour  objet :  1° 
de  verifier  l'exactitude  des  provisions  de  conditions  de  flambage,  calculOes  d’aprOs  la 
mOthode  des  elements  finis,  et  2°  d'Ovaluer  les  avantages  de  l'utilisation  d'un  installation 
d'essai  h  haute  pression  existant  pour  des  essais  sur  module  de  composants  d'appareils 
soumis  &  pression  plus  complexes.  Les  rOsultats  de  1'expOrience  et  ceux  des  calculs  par  la 
mOthode  des  elements  finis  concordaient,  tant  pour  le  mode  que  pour  la  charge  de 
flambage.  Une  autre  mdthode,  visant  k  determiner  la  charge  de  flambage  Olastique  par  la 
mesure  des  frequences  dOcroissantes  des  modes  vibratoires  correspondant  au  modes  de 
flambage,  a  Ogalement  re?u  evaluation  favorable.  Le  reservoir  de  haute  pression  a  bien  jouO 
son  role  dans  l’experience,  bien  que  la  nature  eioignde  de  l'essai  rende  difficile  la  mesure  de 
la  charge  de  flambage. 
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Notation 


shell  radius  to  mid  thickness 


Af  area  of  ring  stiffener 

b  faying  flange  width 

d  distance  from  cylinder  axis  to  frame  centroid 

E  Young’s  modulus 

h  shell  wall  thickness 

L  length  of  shell  bay  between  stiffeners 

n  circumferential  harmonic  wave  number 

PCT  critical  interframe  buckling  pressure 

Pc  pressure  at  which  outer  shell  surface  at  midbay  reaches  yield 

P,  pressure  at  which  inner  shell  surface  at  stiffener  reaches  yield 

gv  yield  stress 

v  Poisson’s  ratio 

pw  density  of  water 

7  A+bh+ 

{[l+»9rt^L]sinh  |  cos  f +11  ~\T^i  '°sh  f  8>n  f  1) 


sinha+sina 


ginha— ain  a 
sinh  a+Bin  a 

cosh  a -cob  a 
Binh  a+sin  a 
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1  Introduction 


The  prediction  of  collapse  loads  for  externally  pressurized  vessels  is  complicated  by  the  un¬ 
certain  effects  of  shape  imperfections,  material  nonlinearities  and  asymmetric  structure  and/or 
load.  Typically,  designs  are  based  on  predictions  of  the  linear  elastic  bifurcation  buckling 
loads  of  ‘perfect’  structure,  together  with  application  of  ‘knockdown  factors’  to  account  for  the 
possible  effects  of  geometric,  material  or  load  uncertainties.  In  lieu  of  the  use  of  knockdown 
factors,  detailed  finite  element  analyses  may,  in  principle,  be  used  to  take  explicit  account  of 
imperfections,  nonlinearities,  and  complex  structure  and  load  functions. 

Detailed  finite  element  analyses  of  submarine  pressure  vessel  structure  are  being  studied  at 
the  Defence  Research  Establishment  Atlantic  (DREA).  This  memo  describes  an  experiment  to 
provide  the  necessary  verification  of  finite  element  interframe  buckling  collapse  predictions  of  a 
ring  stiffened  cylinder  and  to  evaluate  a  high  pressure  tank  facility  for  conducting  model  tests. 

DREA  has  a  high  pressure  water  tank  for  proof  testing  equipment  intended  for  use  in 
deep  ocean  conditions.  To  evaluate  this  facility  for  conducting  pressure  vessel  model  tests,  a 
machined  aluminum  ring  stiffened  cylinder  designed  to  fail  in  an  elastic  interframe  buckling 
mode,  was  tested  to  failure.  The  model  was  machined  from  aluminum  tubing  to  obtain  as 
perfect  a  structure  as  possible  and  then  measured  to  obtain  the  variations  in  shell  thickness 
and  radius.  The  shell  material  was  also  tested  to  obtain  Young’s  modulus  and  yield  stress 
values.  Strain  gauges  were  placed  at  several  locations  on  the  model,  providing  the  only  means 
of  monitoring  its  behaviour  during  the  tests.  The  strain  gauge  information  was  used  to  detect 
the  failure  point  directly  through  observing  sudden  changes  in  strain,  and  indirectly  through 
determining  the  change  in  vibration  frequencies  with  increasing  pressure.  These  test  data  were 
also  compared  to  results  from  an  axisymmetric  finite  element  analysis  carried  out  with  the 
linear  elastic  code  VAST1  developed  by  DREA. 

Section  2  of  this  memorandum  briefly  describes  the  basic  theory  of  buckling  and  collapse 
of  ring  stiffened  pressure  vessels.  The  theory  of  using  vibration  modal  information  *o  nonde¬ 
structive^  determine  buckling  pressure  is  also  discussed.  Section  3  describes  the  test  model 
followed  by  a  discussion  of  the  finite  element  model  of  the  cylinder  in  Section  4.  The  pressure 
tank  facility  and  test  procedure  are  described  in  Section  5  and  the  experimental  and  numerical 
stress  and  collapse  results  are  given  in  Section  6. 


2  Background  Theory 

A  ring  stiffened  cylinder  under  hydrostatic  pressure  (such  as  the  model  investigated  here 
and  shown  in  Figure  1)  can  fail  through  buckling  or  yielding.  Buckling  failure  may  occur  in 
the  elastic  material  phase  of  behaviour  or  after  the  material  has  started  to  yield.  The  critical 
buckling  mode  for  ring  stiffened  cylinders  occurs  in  one  of  two  forms.  One  is  overall  instability 
(Figure  2)  which  usually  occurs  in  a  mode  with  2  to  4  waves  circumferentially  and  a  half  wave 
between  rigid  ends.  The  other  form  is  interframe  instability  (Figure  3)  which  occurs  in  a  mode 
with  a  higher  number  of  circumferential  waves  (6  to  20)  and  half  waves  between  ring  stiffeners. 
In  this  particular  study  interframe  buckling  failure  was  investigated.  The  critical  points  for 
yielding  failure,  shown  in  Figure  4,  are:  the  circumferential  stress  at  the  outer  surface  of  the 
shell  at  midbay  between  stiffeners;  and,  the  axial  stress  at  the  inner  surface  of  the  shell  at  the 
shell-stiffener  connection.  Strain  gauges  were  placed  on  the  model  to  monitor  these  locations. 

The  test  procedure  was,  basically,  to  increa?.,  Uie  pressure  in  the  tank  and  observe  the  strain 
gauge  and  pressure  transducer  responses  for  occurrence  of  buckling  or  yielding. 

Another  method  of  determining  instability  experimentally  is  to  observe  the  values  of  the 
vibration  frequencies  which  correspond  to  buckling  mode  shapes.  For  uniform  models  under 
uniform  loading,  the  mode  shapes  for  vibration  and  buckling  will  be  similar.  It  is  well  known 
that  as  the  stiffness  of  a  structure  is  reduced,  its  natural  frequency  of  vibration  drops  as  the 
square  root  of  the  stiffness  (ie.  frequency  oc  As  the  structure  approaches  its  buckling 

load,  its  stiffness  approaches  zero  and  hence  so  does  its  natural  frequency  of  vibration.  This  is 
easily  understood  for  a  simple  structure  such  as  a  column  under  an  increasing  axial  load  where 
the  principle  buckling  and  vibration  modes  are  well  known.  For  the  ring  stiffened  cylinder,  this 
is  more  complex,  as  there  are  several  buckling  and  vibration  modes  which  may  be  the  principle 
mode  of  buckling  failure.  For  axisymmetric  cylinders,  the  modes  of  buckling  and  vibration  are 
separated  into  circumferential  harmonics,  there  being  an  infinite  number  of  circumferential  har¬ 
monics  each  with  an  infinite  number  of  eigenmodes  describing  the  longitudinal  variation  of  the 
mode  shape.  We  are  concerned  with  only  the  first  eigenmodes  of  each  circumferential  harmonics 
since  they  produce  the  lowest  buckling  load  for  each  harmonic.  Therefore,  in  considering  the 
lowest  buckling  mode  for  each  circumferential  harmonic,  there  is  an  associated  vibration  mode 
whose  frequency  approaches  zero  as  the  buckling  load  for  that  mode  is  reached.  Of  course, 
only  the  frequency  of  one  of  the  harmonics  will  reach  zero,  that  corresponding  to  the  lowest 
(critical)  buckling  mode.  The  frequencies  decrease  as  the  square  root  of  the  applied  load  since 
the  stiffness  terms  for  each  harmonic  are  directly  proportional  to  the  applied  load.  By  using 
measured  vibration  data  from  the  model  as  the  pressure  is  increased,  the  critical  buckling  load 
can  be  determined  by  extrapolating  the  data  to  zero  frequency.  This  technique  has  been  used 
by  Singer2  for  axially  loaded  cylinders.  He  also  proves  that  this  technique  includes  the  effects 
of  imperfections  and  uncertain  boundary  conditions  in  determining  the  actual  buckling  load  of 
the  structure.  In  this  present  study  this  method  is  considered  only  as  a  secondary  investigation, 
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as  the  cylinder  strain  gauge  monitoring  did  not  provide  enough  data  to  distinguish  the  higher 
order  mode  shapes  of  interframe  vibration.  Correlating  the  experimentally  determined  vibra¬ 
tion  frequencies  (from  the  strain  gauge  data)  with  finite  element  predictions  of  the  vibration 
frequencies,  at  increasing  pressure,  allowed  a  means  to  identify  possible  modes,  and  did  give 
some  useful  results. 

Experimental  determination  of  buckling  loads  is  certainly  not  a  new  field,  and  there  have 
been  many  studies  to  investigate  buckling  loads  for  a  variety  of  pressure  vessel  geometries  and 
loads  (reference  3  for  example).  It  has  been  recognized  in  past  studies  that  careful  geometric 
measurement  of  the  model  is  essential  in  interpreting  test  results.  Measurements  were  also 
made  in  this  study,  but  since  the  model  was,  in  fact,  machined  very  near  to  the  ideal  geometry, 
imperfections  were  not  a  significant  factor  in  the  model  failure. 
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3  Test  Model 


The  test  model  was  designed  to  be  as  close  to  a  perfect  structure  as  possible  and  to  fail  in 
an  elastic  interframe  buckling  mode  so  that  the  VAST  linear,  elastic  finite  element  model  would 
represent  its  behaviour.  The  ring  stiffened  model  was  created  by  machining  a  section  of  8x0.5 
inch  aluminum  tubing.  The  dimensions  of  the  model  are  shown  in  Figure  5,  and  a  photograph 
of  the  model  in  Figure  1. 

The  material  was  chosen  from  shop  stock  and  its  exact  type  was  unknown,  although  it  was 
believed  to  be  aluminum  6061-T6.  Uniaxial  tension  and  compression  tests  were  conducted  to 
determine  Young’s  modulus  and  yield  stress.  An  example  of  one  of  the  compression  tests  is 
shown  in  Figure  6.  A  value  of  Young’s  modulus  of  10,000,000  psi,  in  agreement  with  the  usual 
published  values  for  aluminum,  was  derived  from  all  tests.  The  yield  stress  varied  marginally 
from  40,340  to  41,182  psi  at  0.002  strain  over  six  tests.  A  value  of  41,000  psi,  commonly  used 
for  6061-T6,  was  used  as  the  yield  stress  in  this  study. 

The  model  was  manufactured  by  The  Ship  Repair  Unit  (Atlantic)  of  Maritime  Command. 
First,  a  thin  layer  was  machined  off  the  inside  to  ensure  circularity,  and  then  material  was  ma¬ 
chined  off  external  sections  to  create  the  ring  stiffeners.  Dial  gauge  measurements  of  inner  and 
outer  radii  were  made  at  each  ring,  and  at  midbay  between  each  ring,  at  eight  circumferential 
locations.  Figure  7  shows  the  variation  in  the  radii  to  the  outside  of  the  ring  stiffeners  along 
the  length  of  the  shell.  There  are  8  angular  locations,  where  the  symbol  1  is  at  0  degrees.  2  at 
45  degrees,  3  at  90  degrees,  etc.  The  variation  (difference  between  minimum  and  maximum), 
excluding  the  low  value  at  the  cylinder  end,  is  0.038  inches  over  a  nominal  radius  of  4.0  inches. 

Figure  8  shows  the  same  type  of  graph  as  Figure  7  for  the  outside  shell  radii  at  midbay.  The 
variation  is  0.011  inches  over  a  nominal  value  of  mean  (midthickness)  radius  of  3.548  inches  (the 
value  used  in  the  finite  element  model).  Figure  9  shows  the  inside  shell  radii  with  a  variation  of 
0.0068  inches.  Figure  10  shows  the  shell  thickness  with  a  variation  of  0.01  inches  over  a  nominal 
design  value  of  0.06  inches.  From  these  data,  the  maximum  out  of  circularity,  determined  as 
the  maximum  deviation  from  the  mean  radius,  is  0.14  percent  of  the  mean  radius,  which  is  well 
below  recommended  design  values  of  0.5  percent4. 

The  variations  in  the  inner  and  outer  radii  versus  circumferential  angle  were  also  determined. 
Figure  11  shows  this  for  the  center  bay  (where  failure  occured).  A  single  sine  wave  is  evident 
indicating  that  the  measurement  system  may  be  misaligned  or  that  the  model  may  have  a  slight 
overall  bend.  The  thickness  is  given  by  the  difference  between  the  two  lines  which  is  seen  to 
vary  slightly  around  the  model.  In  retrospect,  more  circumferential  locations  should  have  been 
used  to  assess  imperfections  in  higher  circumferential  harmonics. 

Figure  5  also  shows  the  strain  gauge  locations  on  the  inside  of  the  shell  wall.  Both  axial 
and  circumferential  strain  directions  were  measured  at  positions  at  midbay  between  stiffeners 
and  at  the  edge  of  the  stiffener-shell  connection.  The  gauges  were  Micro-Measurement  uniaxial 
and  biaxial  EA-06-125TM-120.  Two  gauges  were  also  placed  at  midbay  in  the  circumferential 


4 


direction  on  the  outside  of  the  shell. 

The  model  was  fitted  with  heavy  endcaps  to  act  as  rigid  ends.  A  solid  aluminum  core,  with 
diameter  one  inch  less  than  the  inner  diameter  of  the  model  was  loosely  fitted  on  the  rigid  ends 
to  act  as  a  stop  mechanism  when  the  model  collapsed.  The  core,  which  was  included  to  avoid 
possible  shock  damage  to  the  pressure  tank  pump  when  collapse  occurred,  did  not  affect  the 
results. 
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4  Analytical  Model 

The  analytical  results  were  determined  by  closed  form  solution  and  numerical  finite  element 
analysis.  The  closed  form  solution  for  the  interframe  buckling  load  is  that  of  Von  Mises  modified 
by  Kendrick5: 


P,r  = 


Eh 


a(n*-  1  +  ffir)1 


+ 


h 2 


[n2(£)2  +  l]2  12a2(J  -  *2 ) 


n 


(1) 


applied  to  a  single  bay  between  frames.  This  formula  has  to  be  used  for  several  values  of  the 
circumferential  harmonic,  n,  to  acheive  the  minimum  value  of  the  buckling  load.  The  closed 
form  solutions  for  the  pressure  values  to  cause  yield  are6or7: 
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m) 


1  +  ~iH 


(2) 


for  the  circumferential  stress  at  the  midbay  outer  fiber  and: 


Ps 


2h°V /  1  X 

a  [1+(*)W 


(3) 


for  the  meridional  (axial)  stress  at  the  inner  surface  at  the  shell-stiffener  connection  (see  Figure 

4). 


The  elastic  interframe  buckling  load  and  stress  values  were  also  determined  by  finite  element 
methods  with  VAST.  Axisymmetric,  three  node  she’.’  elements  were  used  to  model  the  ring 
stiffened  shell  (Figure  12).  One  shell  element  was  used  for  each  stiffener  and  8  elements  were 
used  between  stiffeners.  Hydrostatic  pressure  was  applied  to  the  model  with  an  equivalent 
concentrated  end  load  to  represent  the  pressure  on  the  ends.  The  asymmetric  buckling  loads  and 
vibration  frequencies  are  determined  by  formulating  stiffness  matrices  for  each  circumferential 
harmonic  and  solving  the  eigenvalue  problem  for  each  harmonic8.  The  buckling  eigenvalues  are 
determined  by  solving  for  the  roots  of  the  linear  and  geometric  stiffness  equation: 


[[ K ]  +  \{[Kg]  +  [Kf])}  n  =  0  (4) 

where  A„  are  the  buckling  eigenvalues  directly  proportional  to  the  applied  load,  [ K]n  is  the 
linear  stiffness  matrix  for  circumferential  harmonic,  n,  and  [Kg]„  is  the  geometric  harmonic 
stiffness  matrix.  [Kf]„  is  the  linearized  follower  force  or  load  stiffness  matrix  which  accounts 
for  the  nonlinear  effect  of  the  pressure  vectors  remaining  normal  to  the  shell  surface  during 
deformation.  The  frequency  eigenvalues  are  determined  by  solving  the  similar  problem  for  the 
stiffness  and  mass  matrix: 

[[tf]  +  A[Afjj„  =  0  (5) 
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where  [M]„  is  the  mass  matrix  and  A„  are  the  frequency  eigenvalues.  In  determining  the 
vibration  frequencies  of  the  model  including  the  effects  of  an  applied  hydrostatic  pressure,  the 
eigenvalue  problem  becomes: 

[([*]  + [tfc])  +  A[M]]n  =  0  (6) 

where  [ifcjn  incorporates  the  effect  of  the  applied  prestresss  and  ([if]  +  [ifc?])n  is  the  total 
stiffness  which  reduces  towards  zero  as  the  applied  compressive  pressure  is  increased. 

The  pressure  tank  uses  water  as  its  pressurizing  medsuiTO  so  the  added  fluid  mass  associated 
with  the  cylinder  motion  is  required  in  the  mass  matrix  to  determine  the  vibration  frequencies 
of  the  model  in  the  tank.  At  present  there  are  no  axisyrcmaetric  fluid  elements  in  VAST  which 
can  be  used  in  the  harmonic,  asymmeiiit  mode  shape  determination  (fluid  elements  for  the 
axisymmetric,  n=0,  mode  exist).  Therefore,  the  added  aaasa  rigid  body  motion  of  a  cylinder 
(—  pwiLa2L)  was  determined  and  modeSfed  as  an  eqirwaArrat  density  added  to  the  shell  wall 
material.  This  approximation  should  be  very  good  fen  the  Power  harmonics  (n=0,l,2  and  3), 
but  may  lose  some  accuracy  for  higher  harmonics;  however,  it  was  the  only  method  available 
aside  from  developing  three  dimensional  models,  which  would!  have  been  prohibitively  large  for 
accurate  modelling  of  the  higher  order  harmonic  modes 
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5  Test  Procedure 

The  pressure  tank  is  illustrated  in  Figure  13.  It  has  mainly  been  used  to  proof  test  equipment 
intended  for  use  in  deep  ocean  conditions  to  a  maximum  pressure  of  8500  psi.  The  chamber 
lid  is  fitted  with  a  breach  lock.  Fittings  for  external  monitoring  are  located  in  the  bottom 
and  top  of  the  chamber.  Pressure  is  increased  by  pumping  water  into  the  chamber.  There  is 
no  control  over  the  rate  of  pressure  increase.  Starting  and  stopping  the  pump  was  the  only 
way  to  control  the  tank  pressure.  Since  pressure  increased  fairly  rapidly,  and  since  time  was 
required  to  inspect  the  strain  gauges  and  also  to  collect  enough  data  samples  at  each  load  level 
to  determine  vibration  frequencies,  the  load  was  increased  in  increments  of  50  psi  (by  starting 
and  stopping  the  pump)  until  failure.  The  pump  also  excited  small  shock  waves  in  the  fluid, 
even  while  it  idled.  This  served  as  the  necessary  excitation  function  to  excite  vibration  modes 
in  the  model.  This  characteristic  of  the  apparatus  was  not  known  prior  to  the  experiment,  but 
significant  vibrations  (measurable  from  strain  gauge  data)  occured,  which  spurred  the  secondary 
investigation  into  the  decreasing  frequency  method. 

The  model  was  suspended  vertically  from  a  frame  (Figure  13)  to  minimize  the  strains 
resulting  from  self  weight  and  buoyancy.  Strains  were  measured  in  this  position  and  found  to 
be  negligible.  A  pressure  transducer  with  a  maximum  output  of  1000  psi  was  used  to  measure 
the  applied  pressure. 

Strain  gauge  and  pressure  data  were  recorded  on  analogue  tape.  Chart  recorders  were  used 
to  inspect  the  data  visually  as  the  test  progressed.  A  drop  in  pressure  or  nonlinearities  in  the 
strain  data  were  expected  as  the  model  buckled.  Some  nonlinearity  was  detected  in  a  gauge  on 
the  buckled  bay  but  this  was  only  just  before  the  model  imploded,  which  was  detectable  on  all 
of  the  instrumention,  and  audibly. 

The  vibration  frequencies  were  determined  by  using  a  spectrum  analyzer  on  the  circum¬ 
ferential  strain  gauge  in  the  buckled  bay.  Frequency  spectra  were  generated  for  several  load 
levels. 


8 


6  Results  and  Discussion 


The  test  model  failed  catastrophically  with  no  prior  detectable  indication  of  collapse.  At 
collapse,  the  model  formed  its  buckled  shape,  and  the  shell  completely  tore  away  from  the 
heavy  ring  stiffener.  From  the  stress  and  pressure  data  prior  to  and  during  the  collapse,  shown 
in  Table  1,  the  experimental  collapse  load  appeared  to  be  587  psi.  Figures  14  and  15  show  the 
collapsed  model.  From  Figure  14  it  can  be  seen  that  the  model  failed  in  the  center  bay  (bay 
5,  see  Figure  5)  with  buckling  failure  also  evident  in  bay  4.  Figure  15  shows  the  failed  cross 
section  where  9  circumferential  waves  are  quite  evident.  The  model  collapsed  onto  the  central 
core  in  the  model,  hence  the  flattened  waves.  This  is  quite  clearly  an  interframe  failure  in  mode 
9.  The  actual  elastic  buckling  load  is  probably  a  little  less  than  587  psi  as  the  buckles  should 
have  formed  before  complete  yielding  collapse,  although  this  could  not  be  detected  in  the  strain 
gauge  data.  Prior  to  collapse,  pressure  was  held  steady  at  about  540  psi  for  about  40  seconds 
and  then  increased  towards  590  psi  where  collapse  occured  just  as  the  pump  was  being  turned 
off.  The  pressure  data  in  the  failure  region  are  shown  in  Figure  16.  It  can  be  concluded  from 
this  data  that  the  elastic  buckling  collapse  occurred  somewhere  between  540  and  587  psi. 

The  Von  Mises-Kendricks  formula  (equation  1)  predicted  an  elastic  interframe  buckling 
load  of  556  psi  in  circumferential  harmonic  n=8.  Since  this  formula  assumes  a  section  of  shell 
simply  supported  between  rigid  rings,  it  should  predict  conservative  values  for  a  continuous 
ring  stiffened  cylinder. 

The  VAST  model  shown  in  Figure  12  predicted  a  minimum  buckling  load  of  568  psi  in 
circumferential  harmonic  n=8.  This  was  closely  followed  by  a  buckling  load  of  588  psi  in 
harmonic  n=9.  Figure  17  shows  the  first  two  eigenvalues  for  each  of  the  first  10  circumferential 
harmonics  from  the  VAST  analysis.  This  is  a  typical  graph  for  ring  stiffened  pressure  vessels 
such  as  submarine  hulls,  with  two  minima;  one  for  the  overall  buckling  load  (in  this  case  at 
n=2)  and  one  for  the  interframe  buckling  load  (in  this  case  at  n=8).  Figures  18  and  19  show  the 
VAST  model  mode  shapes  corresponding  to  the  n=8  and  n=9  buckling  harmonics  respectively. 
Note  that  the  mode  for  harmonic  9  (Figure  19)  is  very  similar  to  the  experimental  mode  shape 
(Figures  14  and  15).  This  is  expanded  into  a  three  dimensional  form  in  Figure  20  for  better 
comparison  to  the  experimental  shape.  Since  this  is  a  uniformly  stiffened  cylinder,  it  is  quite 
natural  to  have  several  interframe  modes  with  almost  the  same  buckling  load,  and  the  one  in 
which  it  actually  fails  is  difficult  to  predict.  Hence,  the  fact  that  it  failed  in  harmonic  n=9 
instead  of  the  predicted  n=8  is  neither  unexpected  nor  significant. 

The  strain  gauge  data  were  converted  to  stress  for  some  of  the  gauges.  The  midbay  circum¬ 
ferential  and  meridional  stress  for  the  inner  surface  of  bay  3  and  the  midbay  circumferential 
stress  at  the  inner  and  outer  surface  of  bay  4  are  given  in  Table  2  for  a  section  of  the  test  in  the 
range  of  350  psi  pressure.  Table  1,  which  was  already  briefly  mentioned,  gives  the  meridianal 
stress  on  the  inner  shell  surface  and  the  midbay  circumferential  and  meridianal  stress  on  the 
inner  shell  surface  of  bay  5  (where  failure  occurred)  for  the  test  in  the  failure  range.  Figure 
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21  shows  the  stress  values  for  the  last  80  seconds  of  the  test  corresponding  to  Table  1.  Figure 
16  showed  the  pressure  data  for  the  same  period.  The  influence  of  the  pressure  loading  steps 
can  be  seen  in  the  figures.  As  can  be  seen,  the  data  were  certainly  not  smooth.  Some  of  the 
probable  causes  of  this  were:  60  Hz  line  noise  which  was  very  evident  in  the  frequency  analysis 
of  the  data  (to  be  discussed  later),  pressure  fluctuations  in  the  tank  from  the  pump  vibration, 
and  excitation  of  vibration  modes  in  the  test  model.  It  was  not  expected  that  the  data  would 
fluctuate  so  much;  this  will  have  to  be  investigated  further  in  subsequent  tests. 

In  looking  at  the  measured  midbay  stresses  in  Tables  1  and  2,  it  can  be  seen  that  the 
meridional  values  are  about  half  the  circumferential  values,  as  would  be  expected.  There  is 
very  little  difference  between  the  inner  and  outer  surface  circumferental  stress  indicating  pure 
membrane  behaviour  at  midbay,  and  also  little  difference  between  bays  in  the  central  part  of 
the  model.  These  too,  are  expected  results. 

Agreement  of  the  measured  stress  results  with  the  linear  finite  element  and  formula  predic¬ 
tions  are  not  as  good  as  was  expected.  The  measured  meridional  stress  at  the  stiffener  shell 
connection  at  failure  was  (from  Table  1)  35,699  psi  compared  to  38,350  psi  from  the  finite 
element  results.  The  measured  midbay  stresses  at  failure  were  27,515  psi  in  the  circumferential 
direction  and  13,974  psi  in  the  meridional  direction.  The  corresponding  finite  element  results 
are  33,405  and  16,138  psi,  respectively.  The  Salerno  and  Pulos  formula  for  midbay  circumfer¬ 
ential  stress  gives  34,928  psi.  The  measured  results  seem  to  be  on  the  low  side,  particularly  at 
midbay.  Some  reasons  for  this  may  be: 

•  the  actual  model  is  on  average  slightly  thicker  than  the  0.06  inch  nominal  thickness 

•  the  gauges  may  be  slightly  out  of  line  -  it  would  be  possible  to  calibrate  the  gauges  on 
the  structure  only  by  comparing  them  to  numerical  results,  which  is  the  purpose  of  the 
test,  so  no  calibration  was  done,  and 

•  the  shell  is  experiencing  some  geometric  nonlinearities. 

The  first  two  reasons  should  only  have  small  influences  on  the  stresses.  The  third  reason  is 
possible,  as  the  finite  element  results  indicate  a  radial  deflection  of  0.009  inches  between  midbay 
and  the  ring  stiffener  at  the  failure  pressure.  This  is  0.15  times  the  shell  thickness,  which  is 
within  the  range  where  large  deflection  theory  begins  (0.1  to  0.2  times  the  shell  thickness). 
Large  deflections  (geometric  nonlinearities)  would  increase  the  stiffness  and  thus  lower  the 
strain  gauge  readings. 

The  fact  that  the  meridional  stress  near  the  shell-stiffener  intersection  approaches  the  yield 
stress  of  41,000  psi  explains  the  catastrophic  failure,  and  plentiful  evidence  of  yielding  is  appar¬ 
ent  in  the  collapsed  model.  The  strain  gauges  cannot  detect  the  peak  stress  in  this  area  as  they 
average  over  their  width.  The  finite  element  model  also  averages  over  element  widths.  In  fact, 
a  local  nodal  stress  value,  before  averaging,  of  42,000  psi  occurs  in  the  finite  element  results  at 
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the  connection  point.  It  is  therefore  quite  likely  that  in  the  range  of  590  psi  pressure  a  local 
yield  was  attained  at  the  shell-stiffener  intersection,  leading  to  collapse. 

As  mentioned,  it  was  evident  from  the  strain  gauge  data  that  there  was  sufficient  energy 
in  the  tank  to  excite  the  model  vibration  modes.  This  appears  to  be  a  result  of  unsteady 
motion  of  the  pump  system  producing  shock  waves  in  the  tank.  This  allowed  an  investigation 
using  model  natural  frequency  data  to  predict  buckling  loads.  The  basic  theory  of  using  the 
frequency  squared  versus  the  applied  pressure  has  been  described  in  Section  2.  This  technique 
has  been  investigated  extensively  by  Singer2  for  nondestructive  evaluation  of  buckling  loads 
including  uncertainties  in  boundary  conditions  and  initial  imperfections.  Most  of  his  work 
was  on  axially  loaded,  stringer  stiffened  (stiffeners  in  meridian  axis)  cylinders.  Singer  used  an 
electromagnetic  shaker  to  excite  the  models  and  mode  shapes  were  determined  by  scanning  the 
shell  with  a  microphone.  Vibration  frequencies  for  various  modes  were  measured  over  a  range 
of  applied  loading  up  to  about  60  percent  of  the  predicted  classical  buckling  load.  Straight  lines 
were  extrapolated  through  the  frequency-squared  points  in  the  mode  of  failure  to  the  point  of 
zero  frequency  to  determine  the  buckling  load.  Information  from  numerical  calculations  of  the 
frequencies  versus  applied  load  was  also  used  in  producing  the  curves. 

The  main  drawback  in  applying  this  method  here  was  that  there  was  insufficient  instru¬ 
mentation  to  identify  mode  shapes;  only  frequencies  could  be  determined  from  the  strain  gauge 
data.  Fourier  spectra  for  five  load  levels  were  generated  from  the  circumferential  strain  gauge 
data  in  bay  5,  where  collapse  occurred.  These  spectra  are  shown  in  Figures  22  to  26.  There 
are  some  fairly  strong  peaks  in  addition  to  the  60  Hz  harmonics  in  these  data.  Some  of  the 
peaks  also  shift  with  increasing  applied  load.  This  can  be  more  readily  seen  in  Figure  27  which 
shows  a  cascade  plot  of  frequency  spectra  versus  load  level.  The  predominant  frequencies  were 
squared  and  plotted  against  applied  pressure  in  Figure  28.  Broad  peaks  are  plotted  as  bars 
to  cover  a  range  of  frequencies.  Finite  element  results  for  frequencies  versus  load  were  deter¬ 
mined  and  straight  lines  of  frequency-squared  versus  load  were  plotted  on  Figure  28  for  the 
first  eigenvalues  of  harmonics  0,1, 2, 7, 8, 9  and  10. 

There  are  several  things  to  note  from  Figures  22  to  28.  There  are  frequency  values  of  544  Hz 
and  approximately  470  Hz  that  remain  constant  throughout  the  loading  history.  The  544  Hz 
value  may  correspond  to  a  finite  element  result  of  530  Hz  for  the  first  eigenvalue  of  harmonic  n=0 
which  is  sometimes  referred  to  as  the  ‘breathing’  mode.  This  finite  element  result  remains  at 
530  Hz  throughout  the  increase  in  loading.  This  frequency  shows  a  significant  amount  of  energy 
which  is  another  characteristic  of  the  breathing  mode.  The  470  Hz  value  probably  corresponds 
to  the  finite  element  value  of  418  Hz  for  the  first  eigenvalue  of  the  n=l  harmonic  which  also 
remains  constant  with  increased  load.  This  is  the  fundamental  beam  bending  mode  of  the 
model  which  should  not  change  appreciably  with  increased  pressure.  It  is  also  possible  that 
this  mode  corresponds  to  the  first  eigenvalue  of  the  n=2  harmonic  for  which  the  finite  element 
analyses  give  values  ranging  from  500  Hz  at  50  psi  pressure  to  462  Hz  at  525  psi  pressure.  The 
lines  on  Figure  28  for  harmonics  n=0,l  and  2  pass  through  the  experimental  data  fairly  well. 
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The  experimental  data  points  on  Figure  28  follow  the  general  trend  of  the  numerically 
derived  lines  for  harmonics  n=7,8,9  and  10.  There  are  many  modes  which  exhibit  frequencies 
in  this  range.  Figure  29  shows  the  first  three  eigenvalues  for  each  of  the  first  10  harmonics  at 
a  load  of  50  psi.  Each  point  represents  a  mode  and  there  would  be  more  points  if  more  than 
three  modes  per  harmonic  were  determined.  Therefore,  it  is  difficult  to  ascertain  which  mode 
the  experimental  values  correspond  to.  If  we  assume  that  the  values  correspond  to  the  buckling 
modes,  which  are  the  first  eigenvalues  of  each  harmonic,  identification  is  easier.  This  is  a 
reasonable  assumption  as  these  modes  should  exhibit  the  most  drastic  decrease  with  increased 
load.  Using  Singer’s  method  of  extrapolating  the  data  points  from  frequencies  measured  at 
loads  less  than  60%  of  the  buckling  load,  a  dashed  line  (Figure  28)  is  drawn  through  the  points 
most  likely  to  correspond  to  the  n=9  harmonic  in  which  buckling  failure  occurred.  This  gives  a 
buckling  failure  load  of  555  psi  which  is  within  the  range  of  the  actual  value  which  lies  between 
540  and  587  psi. 

A  combination  of  the  finite  element  analyses  and  limited  modal  data  has  resulted  in  a 
fairly  successful  prediction  of  the  buckling  load.  It  is  planned  to  try  this  again  for  a  model 
which  buckles  in  an  n=2  or  n=3  overall  mode  and  to  identify  the  mode  shapes  with  sufficient 
instrumentation.  If  successful,  this  could  be  useful  in  determining  the  overall  buckling  loads 
of  full  scale  structures,  such  as  submarine  pressure  hulls,  by  measuring  vibrations  at  several 
locations  on  the  circumference  over  the  operating  depth  range  of  the  pressure  hull.  Prior 
determination  of  the  most  likely  buckling  and  frequency  modes  by  a  finite  element  program 
such  as  VAST  would  be  a  necessary  part  of  the  study  to  determine  optimum  positions  for 
instrumentation  and  to  evaluate  the  results. 

Another  aspect  of  vibration  and  stability  analysis  is  that  of  parametric  resonance.  This 
occurs  when  a  periodic  forcing  function  drives  certain  vibration  frequencies.  The  Mathieu9 
equation  describes  this  behaviour  where  the  loading  function  is  a  coefficient  of  the  displacement 
in  the  equation  of  motion.  For  certain  driving  frequencies  and  amplitudes  of  the  forcing  function, 
the  resultant  vibration  motion  can  be  significantly  amplified  to  the  point  of  instability  in  the 
corresponding  buckling  mode.  Periodic  loading  of  amplitudes  less  than  the  static  collapse  load 
can  cause  failure  in  this  manner.  It  is  possible  in  this  case  that  periodic  loading  from  the  pump 
system  may  have  excited  the  n=9  mode  near  the  failure  point  and  caused  the  final  pressure 
perturbation  necessary  for  collapse. 
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7  Concluding  Remarks 

The  experimental  determination  of  the  interframe  buckling  load  of  a  ring  stiffened  cylinder 
served  several  purposes.  The  finite  element  code  VAST  successfully  predicted  the  buckling  load 
and  mode,  and  also  some  of  the  predominant  vibration  frequencies.  The  fact  that  the  predicted 
buckling  load  for  the  perfect  model  was  very  close  to  the  actual  buckling  load  is  largely  due 
to  the  excellent  machining  of  the  model  by  the  Ship  Repair  Unit  (Atlantic)  which  resulted  in 
minimal  imperfections  in  the  model.  More  extensive  measurement  of  the  initial  geometry  would 
have  been  desirable  to  detect  imperfections  in  the  higher  order  modes  where  buckling  occurred. 

The  pressure  tank  facility  was  shown  to  be  suitable  for  model  buckling  tests  although  the 
inability  to  see  the  model  during  the  test  or  monitor  general  displacements  in  some  manner 
made  determination  of  the  exact  buckling  load  difficult.  The  comparison  of  measured  and 
predicted  strain  data  was  not  as  good  as  might  be  expected,  but  was  still  reasonable.  The 
measured  data  showed  the  expected  stress  distribution. 

The  investigation  of  the  decreasing  frequency  method  for  nondestructively  determining 
buckling  loads  showed  considerable  promise,  and  will  benefit  from  further  work  in  this  di¬ 
rection.  More  effort  in  measuring  vibration  mode  shapes  as  opposed  to  just  the  frequency 
values  should  make  this  a  viable  technique  to  employ  on  full  scale  structures  such  as  opera¬ 
tional  submarines.  Determination  of  modes  in  shell  structures  which  have  many  modes  in  the 
same  frequency  range  will  prove  difficult,  and  it  is  unlikely  that  the  interframe  modes  of  a  full 
scale  structure  could  be  definitively  measured;  however,  overall  buckling  and  vibration  modes 
with  circumferential  harmonics  of  n=2  or  n=3  may  be  discernable. 
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TIME  SIGE8  SIGC5  SIGM5  PRESSURE 

(SEC)  (PSI)  (PSD  (PSD  (PSD 

75.8180  -32294.08  -26382.28  -13208.72  545.8 

75.9850  -31525.17  -25613.38  -12954.98  543.3 

76.1520  -32403.92  -26378.99  -14088.57  536.0  SIGE8 : 

76.3190  -31635.01  -25869.31  -13137.32  545.8  MERIDIANAL  STRESS 

76.4860  -31635.01  -25942.91  -13063.72  543.3  AT  SHELL- STIFFENER 

76.8530  -32294.08  -25721.02  -13577.79  536.0  CONNECTION  IN  BAY  8 

76.8200  -31415.33  -26236.19  -13062.62  543.3 

76.9870  -31744.86  -25759.47  -13101.07  545.8  SIGC5 : 

77.1540  -32184.23  -26013.21  -13869.98  528.7  CIRCUMFERENTIAL 

77.3210  -31635.01  -26126.35  -13026.38  545.8  STRESS  IN  BAY  5 

77.4880  -31964.55  -25685.87  -13174.67  543.3 

77.6550  -32623.61  -27040.25  -13719.49  545.8  SIGM5 : 

77.8220  -32513.76  -28649.46  -14739.94  545.8  MERIDIAN AL 

77.9690  -32403.92  -28466.02  -14777.29  548.2  STRESS  IN  BAY  5 

78.1560  -32184.23  -26381.19  -13502.00  550.6 

78.3230  -31964.55  -25906.66  -12953.88  545.8 

78.4900  -32294.08  -26125.25  -13319.66  543.3 

78.6570  -32823.61  -26271.34  -13465.75  550.6 

78.8240  -31964.66  -25942.91  -13063.72  553.1 

78.9910  -32403.92  -26271.34  -13465.75  545.8 

79.1580  -33172.83  -26599.77  -13867.78  548.2 

79.3250  -32623.61  -26748.06  -13427.31  560.4 

79.4920  -32513.76  -26454.78  -13428.40  558.0 

79.6590  -32623.61  -26637.12  -13684.34  545.8 

79.8260  -33172.83  -26820.56  -13646.99  562.8 

79.9930  -32733.45  -26308.69  -13282.31  565.3 

80.1600  -32953.14  -26637.12  -13634.34  548.2 

80.3270  -33612.20  -27404.93  -14231.36  562.8 

80.4940  -32733.45  -26748.06  -13427.31  567.7 

80.6610  -33062.98  -26857.91  -13463.55  565.3 

80.8280  -33172.83  -26784.31  -13537.15  565.3 

80.9950  -33502.36  -26673.37  -13794.18  555.5 

81.1620  -33722.05  -27111.65  -14232.46  558.0 

81.3290  -34490.95  -27221.49  -14268.71  567.7 

81.4960  -33941.73  -27258.84  -14085.27  570.2 

81.6630  -34490.95  -27442.28  -14047.92  579.9 

'  81.8300  -34490.95  -27368.68  -14121.52  575.0 

81.9970  -34600.80  -27185.24  -14158.87  577.5 

82.1640  -34600.80  -27185.24  -14158.87  575.0 

82.3310  -34381.11  -26929.31  -13976.53  560.4 

82.4980  -34490.95  -26895.25  -13280.12  560.4 

82.6650  -34710.64  -27004.00  -13609.65  579.9 

82.8320  -34271.27  -27187.44  -13572.30  579.9 

82.9990  -34381.11  -26895.25  -13280.12  587.2 

83.1660  -34710.64  -26821.66  -13353.71  589.7 

83.3330  -34820.48  -27590.57  -13607.45  587 . 27FAILURE 

63.5000  -35699.23  -27515.87  -13974.33  589. 7) 

83.6670  -35259.86  -26637.12  -13684.34  579.9 

83.8340  -15817.51  -23255.03  -16189.88  506.7 

84.0010  111601.30  21225.12  -97631.37  6.5 

84.1680  98639.73  43713.44  -92069.98  50.4 


TABLE  1:  Strain  and  Pressure  Data  Near  Failure  Point 
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TIKE 

(SEC) 

0.0000 

0.1670 

0.3340 

0.5010 

0.6680 

0.8350 

1.0020 

1.1690 

1.3360 

1.5030 

1.6700 

1.8370 

2.0040 

2.1710 

2.3380 

2.5050 

2.6720 

2.8390 

3.0060 

3.1730 

3.3400 

3.5070 

3.6740 

3.8410 

4.0080 

4.1750 

4.3420 

4.5090 

4.6760 

4.8430 

5.0100 

5.1770 

5.3440 

5.5110 

5.6780 

5.8450 

6.0120 

6.1790 

6.3460 

6.5130 

6.6800 

6.8470 

7.0140 

7.1810 

7.3480 
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TABLE  2:  Strain  and  Pressure  Data  Near  350  psi 
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Figure  1:  Photograph  of  Model  Ring  Stiffened  Pressure  Vessel 
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Figure  2:  Overall  Buckling  Mode  for  A  Ring  Stiffened  Cylinder 
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Figure  3:  Interframe  Buckling  Mode  for  A  Ring  Stiffened  Cylinder 
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Figure  4:  Critical  Stress  Locations  for  A  Ring  Stiffened  Pressure  Vessel 
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Figure  5:  Schematic  of  Model  Ring  Stiffened  Cylinder 
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Figure  6:  Measured  Stress-Strain  Curve  for  Model  Material 


Figure  7:  Radii  to  Outer  Ring  Surfaces  over  Length  and  Circumference 
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Figure  8:  Radii  to  Outer  Shell  Surface  over  Length  and  Circumferenc* 


Figure  9:  Radii  to  Inner  Shell  Surface  over  Length  and  Circumference 
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Figure  10:  Shell  Thickness  over  Length  and  Circumference 


Figure  11:  Radii  of  Shell  Wall  for  Axial  Location  13  Inches  over  Circumference 
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Figure  12:  Axisymmetric  Finite  Element  Model  of  Ring  Stiffened  Pressure  Vessel 
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Figure  13:  Schematic  of  Pressure  Tank  Facility  and  Test  Setup 
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Figure  17:  Buckling  Load  versus  Circumferential  Harmonic  from  Finite  Element  Results 
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Figure  18:  Buckling  Mode  in  Harmonic  n=8  from  Finite  Element  Analysis,  PeT— 568  psi 
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Figure  19:  Buckling  Mode  in  Harmonic  n=9  from  Finite  Element  Analysis,  Pcr=588  psi 
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Figure  20:  Three  Dimensional  Expansion  of  n=9  Buckling  Mode  from  Finite  Element  Analysis 


Figure  21:  Stress  Data  in  Region  of  Failure 
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Figure  22:  Fourier  Spectra  from  Gauge  in  Failed  Bay  at  50  psi  Pressure 
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Figure  23:  Fourier  Spectra  from  Gauge  in  Failed  Bay  at  275  psi  Pressure 


Figure  24:  Fourier  Spectra  from  Gauge  in  Failed  Bay  at  325  psi  Pressure 
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Figure  25:  Fourier  Spectra  from  Gauge  in  Failed  Bay  at  475  psi  Pressure 
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Figure  26:  Fourier  Spectra  from  Gauge  in  Failed  Bay  at  525  psi  Pressure 
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Figure  27:  Cascade  Fourier  Spectra  from  Gauge  in  Failed  Bay  for  Pressures  Ranging  from  50 
psi  to  525  psi 
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Figure  28:  Frequency-Squared  versus  Applied  Pressure  from  Gauge  in  Failed  Bay  and  Results 
of  Finite  Element  Analyses 
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Figure  29:  Natural  Frequency  versus  Circumferential  Harmonic  from  Finite  Element  Results 
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